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ABSTRACT 


First,  the  analytical  technique  is  shown  to  be  vaiid  for  the  computation  of  the  statistics 
of  a  blade's  vibratory  amplitude  when  the  distributions  of  modal  parameters  of  a  mistuned 
bladed  disk  assembly  are  non-Gaussian.  The  results  from  the  analytical  technique  are 
compared  with  those  from  numerical  simulations  for  triangular  and  uniform,  distributions.  It 
has  been  found  that  the  probability  density  function  of  tne  amplitude  is  insensitive  to  the  types 
fo  mistuning  distributions. 

Second,  an  analytical  technique  has  been  developed  to  efficiently  compute  the  probability 
density  function  of  tne  maximum  ampitiude  on  a  mistuned  bladed  disk  assembly.  This  technique 
uses  the  direct  Taylor  series  expansion  in  terms  of  the  perturbation  in  an  amDlitude  as  a 
function  of  perturbations  in  modal  stiffnesses.  The  validity  of  the  technique  has  been 
corroborated  by  comparison  with  the  results  from  numerical  simulations. 

Lastly,  the  statistics  of  the  forced  response  of  a  structurally  and  aerodynamically 
coupled  bladed  disk  assembly  have  been  computed  efficiently  by  the  analytical  technique.  The 
results  from  the  analytical  technique  agree  well  with  those  from  numerical  simulations.  The 
effects  of  the  following  parametrs  on  the  statistics  of  the  maximum  amplitude  have  been 
studied:  the  aerodynamic  couplings  among  blades,  the  fluid  density  and  the  cascade  stagger  angle. 
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SECTION  1 


PROBABILISTIC  ANALYSIS  OF  FORCED  RESPONSE  OF  A  PLADED 
DISK  ASSEMBLY  WITH  VARIOUS  MISTUNING  DISTRIBUTIONS 
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ABSTRACT 

Ir,  this  paper,  efficient  algorithms  have  been  developed  to  compute  the 
statistics  of  forced  response  of  a  mistuned  bladed  disk  assembly.  The 
technique  is  valid  for  Gaussian  and  non-Gaussian  distributions  of  modal 
parameters  as  well.  The  validity  of  this  technique  is  shown  by  comparison 
with  results  from  numerical  simulations  for  triangular  and  uniform 
distributions.  The  results  for  Gaussian,  triangular  and  uniform 
distributions  indicate  that  the  cumulative  probability  density  function  of 
a  blade's  amplitude  primarily  depends  on  the  standard  deviations  of  modal 
parameters  rather  than  on  the  type  of  distributions. 


NOMENCLATURE 


As 

o 

aij 

B 

c 

E(  ) 
f 

91  (  ) ,  92(  ) 
H 

He 

i 

I 

h 

mt 


amplituue  of  the  jth  mass  (mistuned  systei) 

elements  of  matrix  H~^ 

amplitude  of  each  blade  in  tuned  system 

damping  coefficient 

expected  value 

external  force  vector 

probability  density  function 

matrix  with  elements  composed  of  m^,  k^,  c  and  u> 

complex  matrix  corresponding  to  H 

/  -1 

identity  matrix 

modal  stiffness  of  each  blade  (tuned  system) 
modal  mass  of  each  blade  (tuned  system) 
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N  number  of  blades 

NS  number  of  simulations 

r  order  of  engine  excitation 


R  A/kt 

U  vector  defined  by  Eq.  (10) 

uj  jth  element  of  U 

V  vector  defined  by  Eq.  (11) 

Yj  jth  element  of  V 

X  response  vector  (mistuned  system) 

Xt  resnonse  vector  (tuned  system) 

yj  correlation  coefficient  between  U2j  and  U2j-i 

£  damping  factor  (-  c/(2mt  m)  for  resonance  condition) 

6Aj  amplitude  of  <5x j 

6X  X-Xt 


6xj 

jth  element  of  <5X 

phase  of  <5xj 

A 

standard  deviation  of  <Skj 

standard  deviation  of  uj 

phase  of  Xj 

a] 

excitation  frequency 

1.  Introduction 

Mistuning  refers  to  variation  in  modal  properties  of  blades  which 


is 


caused  by  minor  differences  in  their  geometries  that  arise  during  the 
manufacturing  process.  It  has  received  wide  attention  in  the  existing 
literature  [1,2]  because  even  a  small  amount  of  mistuning  can  lead  to  a 
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large  variation  in  the  blade's  amplitudes  within  the  same  assembly.  In 
fact,  the  amplitude  of  vibration  of  the  worst  blade  may  often  be  as  high  as 
two  times  the  amplitude  of  a  blade  of  a  perfectly  tuned  system.  Since  the 
deviations  of  modal  properties  of  blades  from  their  nominal  values  are 
randomly  distributed,  it  is  important  for  the  designer  to  predict  the 
statistics  of  forced  response,  particularly  the  probability  that  the 
maximum  amplitude  on  any  disk  would  exceed  a  critical  value,  in  a 
computationally  efficient  manner.  For  this  purpose,  an  efficient 
computational  technique  has  been  presented  by  the  authors  in  their  earlier 
papers  [3-4]  where  the  distributions  of  modal  parameters  are  assumed  to  be 
Gaussian.  In  this  paper,  the  analytical  technique  is  extended  to  include 
non-Gaussian  distributions.  Also,  an  efficient  algorithm  has  been 
developed  to  invert  the  tuned  system  matrix  which  further  reduces  the 
computation  time. 

The  distributions  of  modal  parameters  are  often  assumed  to  be  Gaussian 
for  the  sake  of  analytical  simplicity.  However,  in  practice,  they  may  net 
be  Gaussian.  Therefore,  the  following  issue  arises:  does  one  type  of 
distributions  of  modal  parameters  leads  to  worse  scatter  in  amplitudes  than 
another?  This  issue  has  been  addressed  in  this  paper  by  studying  the 
distributions  of  amplitudes  for  three  types  of  distributions  of  modal 
parameters  with  the  same  standard  deviations:  Gaussian,  uniform  and 
triangular.  It  is  found  that  the  cumulative  probability  distribution 
functions  of  a  blade's  amplitude  is  almost  identical  for  these  three  types 
of  parameters'  distributions. 
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The  model  of  bladed  disk  assembly  is  chosen  to  be  as  shown  in  Figure  1. 
Even  though  it  is  a  simple  model,  it  contains  all  the  basic  characteristics 
of  a  bladed  disk  assembly  and  has  been  used  by  the  authors  in  their 
previous  papers  [3-5].  In  Figure  1,  mt  and  kj  are  the  modal  mass  and 
stiffness  of  each  blade,  respectively.  The  phenomenon  of  mistuning  is 
simulated  by  considering  kj  as  a  random  variable.  The  stiffness  kC 
represents  the  coupling  between  adjacent  blades  due  to  disk's  flexibility. 

First,  the  technique  for  the  efficient  computation  of  the  statistics 
cf  forced  response  is  described.  The  accuracy  of  this  technique  is  then 
studied  for  triangular  and  uniform  distributions  by  comparison  with  results 
from  numerical  simulations.  Lastly,  the  values  of  probability  that  a 
blade's  amplitude  exceeds  certain  critical  value  are  examined  as  a  function 
of  excitation  frequency  for  Gaussian,  triangular  and  uniform  distributions 
of  k  j  . 

2.  Computation  of  Statistics  of  Forced  Response 

The  computational  technique  is  presented  for  the  system  shown  in 
Figure  1.  However,  it  can  be  directly  applied  to  a  more  complex 
spring-mass  model  of  a  bladed  disk  assembly;  e.g.,  the  model  developed  by 
Griffin  and  Hoosac  [6]. 

The  system  of  differential  equations  governing  the  dynamics  of  a 
mistuned  bladed  disk  model  depicted  in  Figure  1  is  described  as  follows: 


m,  x.  +  k.x.  +  kC(x.-x.  ,)  +  kC(x.-x.  ,)  +  cx. 
t  J  J  J  {  J  j+1'  v  J  j-1'  J 


f0  cos(u>t+i|>.) ;  j-1,  2 


N 


(1) 
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where  N  is  the  number  of  blades.  The  excitation  phase  ijjj  is  defined  for  a 
rth  engine  order  excitation  as  follows: 

-  2n  r(j-l)/N;  r  -  0,  1,  2,  ....  N-l  (2) 

Represent 

kj  -  kt  ♦  6k j  (3) 

where  k^  is  the  stiffness  of  the  perfectly  tuned  system.  The  deviation  6kj 
is  taken  to  be  a  ranacm  variable  with  zero  mean;  i.e.,  E(6kj)  -  0.  Also, 
the  distribution  of  6kp  is  assumed  to  be  independent  of  6kj(p*j).  It 
should  be  noted  that  there  is  no  assumption  made  on  the  nature  of  the 
distribution  of  6k j . 

Now  the  response  (x^j)  of  the  tuned  system  satisfies  the  following 
equation: 


Vtj  *  ktxtj  *  kc<Vxt(M)>  *  kc(xtj-xt(j-i)> 

+  c  x,.  -  f  cos(ut+i|».)  (4) 

tj  0  J 

Subtracting  (4)  from  (1), 

m.  6x  .  +  (k.  +  6k  .)  6x  .  +  kC(6x  .-5x  .  .)  +  kC(6x  .-6x  .  .) 

^  J  ^  J  J  J  J +  ^ 

*  c  6xj  ■  -{kj  xtj  (5) 

where  6xj  -  xj  -  x^j. 

In  steady  state,  the  response  xj  and  x^j  are  sinusoidal.  Therefore,  6xj 
would  also  be  sinusoidal.  Let 


and 


xj  -  Aj  cos  (uit+<f>j) 
x^j  -  B  cos  (ait+Qj ) 
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(6) 

(7) 


6x  -j  -  6Aj  cos  (ut+6<j)j)  (8) 

Substituting  (7)  and  (8)  into  (5)  and  then  equating  the  coefficients  of 
cos  ait  and  sin  tut  on  both  sides,  the  following  system  of  2N  algebraic 
equations  is  obtained: 


(H+<5H)U  -  -5HV 

(9) 

where 

U  -  [6Ai  ccs6<j>i 

5Aj  si n64>i  •  •  •  •  dAf^  cos6^ 

6Am  sin64>N]T 

(10) 

V  -  [ B  cosGi 

B  sinGi  .  .  .  .  B  cos0fj 

B  sin0fj]T 

(11) 

The  matrices  H  and  6H  are  defined  in  the  Appendix  A.  Since  the  matrix  H  is 
non-singular  for  a  damped  system,  the  equation  (9)  yields 
U  -  -( I+H-16H)-1  H-^HV 

-  -(H~l6H  -  (HdH)2  +  (H6H)3  -  .  .  .}V  (12) 

Therefore,  the  expression  for  each  element  of  the  vector  U  has  the 
following  form. 


N  N  N  N  N 

*,  +  l  S  w  j  6k. 6k  -  l  l  l  J  |  ,  6k ,6k  6k 

J  j-l  s-1  PJ,S  J  s  j-1  s-1  i-1  J  s  x 


+  .  .  .  .  (13) 

where 

tp.j  '  aP,2j-l  v 2 j — 1  +  ap,2j  v2j  (14) 

wP,j,s  "  ap »2j— 1  t2j-l,s  +  ap ,2j  t2j,s  (15) 

*P,j,s,i-  aP,2j-l  w2j— 1  ,s , £  +  ap ,2j  w2j  ,s , i  (lfi) 

aP»j  *  ^ H_1)p , j 
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An  efficient  algorithm  has  been  outlined  in  Appendix  B  for  the  inversion  of 
matrix  H.  Therefore,  all  the  coefficients  in  (13)  can  be  computed  in  an 
efficient  manner. 

The  equation  (13)  indicates  that  each  Up  is  a  linear  combination  of  a 
large  number  of  random  variables  which  are  not  highly  dependent  [4], 
Therefore,  the  central  limit  theorem  [7]  indicates  that  the  distribution  of 
Up  can  be  approximated  to  be  Gaussian. 

The  number  of  terms  in  the  series  expansion  (12)  which  make 
significant  contributions  depends  on  the  value  of  R/£  where  R  -  /E(6kj2)/kt 
and  5  is  the  damping  facto*'.  When  R/J,  is  small,  the  higher-order  terms  in 
(13)  can  be  neglected  and  Up  can  he  expressed  a: 


u 


P 


N 

-.1 

j"1 


(18) 


When  the  distribution  of  <5kj  is  Gaussian,  the  distribution  of  Up  is  exactly 
Gaussian  for  any  number  of  blades  [3].  However,  for  non-Gaussian  6k j ,  the 
distribution  of  Up  can  be  approximated  to  be  Gaussian  only  for  a  large 
value  of  N. 

Using  central  limit  theorem,  the  distribution  of  e  (-  n  u 2 j - 1  +  U  u 2 j ) 
can  also  be  approximated  to  be  Gaussian  for  any  q  and  y.  Therefore,  U£j_i 
and  U2j  can  be  described  as  jointly  Gaussian  [8],  i.e., 


9 1 ( u 2 j— 1  *  u2j )  “ 


2^-1^  (1"V 


Tjn 


exp  {- 


1 


ZO-Yp  °2j-l 


2y.ti3 


t-i-V])  (19) 


8 


where 

a  -  u2j_i  -  E(u2j-i) 

8  -  u2j  -  E(u2j } 

ar2  -  E(up  -  (E(ur))2 

Tj  “  (E^U2j-lU2j^  ~  E^u2j-1)  E^u2j  ^ ) 7 C °2j-l  ’  a2j5 

From  (19),  the  probability  density  function  of  a  blade's  amplitude  can  be 
calculated  as 


where 


92(Aj)  - 


2„„2.-ia2j.  /(1-Yj) 


2tt 

exp  [- 
o 


_I1»'  1  dv 


(20) 


2 

(A.cosv-aJ 
Z(v)  .  Xg  — i_ 

Vl 


2yi(A.cosv-a1)  (A.sinv-M 

—  J.,  V  .  *  - - - 


(A -sinv-0, ) 2 
+  — 

a2j 

°1  “  V 2 j  —  1  +  E^2j-l5 
Bl  -v2.  +  E(u2j) 
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The  probability,  Pr,  that  Aj  will  be  less  than  a  critical  value,  Ac,  can  be 
calculated  using  (20)  as  follows; 


It  should  be  noted  that  the  mean  and  variance  of  Up  can  be  calculated 
without  performing  numerical  simulations.  However,  the  number  of  terms  in 
(13)  which  would  yield  accurate  results  has  to  be  determined.  In  this 
paper,  E(up)  and  E(up2)  are  computed  (Appendix  III)  on  the  basis  of  first 
three  terms  in  (13)  for  three  types  of  distributions  for  <5kj:  Gaussian 
[4],  uniform  and  triangular.  It  has  been  found  that  the  results  are  fairly 
accurate  at  a  realistic  value  of  R/£. 

3.  Results  and  Discussions 

3.1  Accuracy  of  the  Analytical  Technique  for  Non-Gaussian  Distributions 
In  order  to  verify  the  accuracy  of  the  analytical  technique  for 
non-Gaussian  distributions  of  modal  stiffness  5kj,  the  results  from 
numerical  simulations  are  compared  to  those  obtained  analytically.  Two 
types  of  distributions,  triangular  (Figure  ?.)  and  uniform.  (Figure  3),  are 
considered  primarily  because  they  can  be  conveniently  generated  using  IMSL 
subroutines  [9]  GGTRA  and  GGU8S  for  numerical  simulations.  These  IMSL 
subroutines  are  used  to  choose  the  modal  stiffness  of  each  blade  from 
populations  with  specified  mean  and  variance.  The  response  of  the 
resulting  bladed  disk  assembly  is  obtained  by  solving  the  linear  system  of 
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algebraic  equations  [3]  in  terms  of  Aj  cos <J)j  and  Aj  s i n j .  The  statistics 
of  A-j  cos 4>j ,  Aj  sin 4> j  and  Aj  are  then  generated  by  considering  a  large 
number  (arbitrarily  chosen  to  be  500)  of  bladed  disk  assemblies  with 
different  sets  of  blades'  modal  stiffnesses. 

The  nominal  values  of  modal  parameters  are  chosen  from  the  paper  by 
Griffin  and  Sinha  [5]  and  are  presented  in  Table  1.  The  number  of  blades 
is  selected  to  be  twenty  four.  The  excitation  frequency  w  is  chosen  to 
represent  resonance  condition.  The  standard  deviation  (.\)  of  each  blade  is 
taken  to  be  8,00C.  At  this  value  of  A,  the  ratio  of  standard  deviation  and 
mean  value  of  individual  blade's  natural  frequency  is  found  to  be  0.0096 
which  represents  the  typical  value  found  in  practice.  The  damping 
coefficient  c  is  chosen  to  be  1.443.  This  value  of  damping  corresponds  to 
about  1%  damping  factor  which  is  typically  observed  in  engines  from 
aerodynamic  and  structural  sources. 

The  results  for  triangular  distribution  are  presented  in  Figures  4-6 
for  third  engine  order  excitation.  The  distributions  of  Aj  cos<f>-j  and 
Aj  sin<j>j,  obtained  from  numerical  simulations,  are  examined  for  each  blade. 
Each  of  these  distributions  is  compared  to  the  Gaussian  distribution  with 
same  mean  and  variance;  e.g..  Figures  4  and  5.  The  chi-square  test  of 
hypothesis  [10]  indicates  that  the  distributions  of  Aj  cos 4>j  and  Aj  sin<j>j 
are  Gaussian,  which  confirms  the  applicability  of  central  limit  theorem. 

The  distributions  of  amplitudes,  obtained  "from  eq.  (20),  are  then  compared 
to  those  from  numerical  simulations;  e.g.,  Figure  6.  The  chi-square  test 
of  hypothesis  again  indicates  that  the  probability  density  function  of 
amplitude  predicted  by  the  analytical  technique  is  accurate. 
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Next  the  results  for  uniform  distribution  are  presented  in  Figure  7-9, 
Again,  the  distributions  of  Aj  costfij  and  Aj  sin^j  are  found  to  be  Gaussian, 
Figures  7  and  8.  Also,  the  probability  distribution  functions  of 
amplitudes,  obtained  from  eq.  (20),  compare  well  with  tKose  from  numerical 
simulations;  e.g.,  Figure  9. 

3.2  Probability  That  The  Amplitude  is  Greater  Than  A  Critical  Value 

Using  the  equation  (21),  the  probability  (Pr)  that  a  blade's  amplitude 
exceeds  a  critical  value  can  be  easily  calculated.  Assuming  that  the 
critical  value  of  the  amplitude  is  1.4 B,  the, values  of  Pr  are  shown  in 
Figure  10  as  a  function  of  excitation  frequencies  for  all  the  three  types 
of  distributions  of  6k j  and  three  different  engine  order  excitations,  it 
is  found  that  the  maximum  value  of  Pr  for  each  engine  order  excitation 
occurs  at  a  frequency  which  is  slightly  below  the  corresponding  undamped 
natural  frequency.  Also,  the  values  of  Pr  are  almost  identical  for  three 
different  distributions  of  6k j .  These  results  indicate  that  the 
probability  that  the  amplitude  would  exceed  a  certain  value  depends 
primarily  on  the  standard  deviation  of  6k j ;  e.g.  Figure  11. 

4.  Conclusions 

In  this  paper,  the  analytical  technique  developed  in  [4]  is  shown  to 
be  valid  for  non-Gaussian  distributions  of  modal  parameters.  This 
technique  provides  the  probability  distribution  function  of  each  blade; 
i.e.,  the  probability  that  a  blade's  amplitude  is  greater  than  a  certain 
critical  value  can  be  easily  calculated.  The  validity  of  this  technique  is 
established  by  comparison  with  results  from  numerical  simulations  for 
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triangular  and  uniform  distributions  of  blades'  medal  stiffnesses.  Even 
though  the  analytical  technique  is  presented  for  a  simple  bladed  disk 
assembly  (Figure  1),  it  can  be  directly  applied  to  a  more  complex 
spring-mass  model  of  a  bladed  disk  assembly;  e.g.,  the  model  used  in  [6]. 
Also,  an  efficient  algorithm  which  uses  fast  Fourier  transform  has  been 
developed  to  invert  the  tuned  system  matrix  which  further  reduces  the 
computation  time. 

The  values  of  probability  that  an  amplitude  is  greater  than  a  certain 
critical  value  are  found  to  be  almost  identical  for  Gaussian,  uniform  and 
triangular  distributions  of  modal  stiffness  with  the  same  values  for 
standard  deviations.  In  other  words,  the  cumulative  probability  density 
function  of  a  blade's  amplitude  depends  primarily  on  the  standard 
deviation  of  modal  stiffness. 
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Appendix  A 
Matrices  H  and  6H 
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hi  h3 
-h3  hi 

-h2  0 
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0  -h2  -h2  G 
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hl  h3 
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hl 


kt  +  2kC  -  mto)2,  h2 


kC,  h3 


-CU) 


6k  i 


6k  i 


6k2 


6k  2 


6kN 


(A— 2) 


6k  N 
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Appendix  B 

Inversion  of  Matrix  H 

It  can  be  seen  from  eq.  (A-l)  that  the  matrix  H  corresponds  to  a  perfectly 
tuned  system  and  is  circular  in  nature.  In  order  to  obtain  an  efficient 
technique  to  determine  H"*,  it  is  convenient  to  use  complex  excitation;  i.e., 
the  differential  equation  for  the  tuned  system  is  written  as  follows 

mt  xy  *  kt  xtj  *  kc(xtj  -  xt(j*D>  ♦  kC<xu  -  xt{j-u) 


(8-1) 


The  steady-state  response  x^j  can  be  represented  as  follows: 


xt.  -Ajeltut  (B-2) 

Substituting  (B-2)  into  (B-l)  and  equating  the  coefficients  of  e’ut  on  both 
sides,  the  following  system  of  N  algebraic  equations  is  obtained: 


Hc5?  ’  -r 


( 8-3) 


where 


( B-4) 


c 


-kC 


0 


kk 


kk  -  (kfc  +  2kC  -  +  icw 

A*r  -  [A*ifr  A*2,r  •  • 
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•  •  •  A*N,r3^ 


(B-5) 


i>,  i\p2 
f*  «  f„[e  1  e  L 

-T  C 


.  .e'V 


( B-6) 


For  the  rth  engine-order  excitation, 


ijjj  -  2nr(j-l)/N 


(B-7) 


From  ( B-3) , 


A*  -  H/1  f 
-r  c  -r 


( B— 8) 


Representing 


(H  _1)  .  -  a*  . 

'  c  'p,j  p,j 


and  using  (B-6),  (B-7)  and  (B-8), 


(B-9) 


A*1  r  "  -  a*l  '  el2irr(i~1]/N  ;  r  -  0,  1,  2,  .  .  .,  N-l  ( B-10) 
’  i*l  ’ 

The  equation  (B-10)  indicates  that  A*i>r  'is  discrete  Four:,er  transform  [11]  of 
a*l,t'  Therefo,  ;„k 


M  "  R 


A*  — 12rrr  ( £— 1)  /N.  ,  ,  - 


N 


(  B-ll) 
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For  each  engine-order  excitation,  A *i>r  can  be  computed  [3]  by  analyzing  an 
equivalent  single-degree  of  freedom  system  and  is  given  by 


1,r  (ke  -  cV 


ctuf 


(K 


TZ 


— 2 — T~ 
to  m^}  +  (.‘'to 


(B-12) 


where 


kg  -  kt  +  4kC  sin2(jrr/N) . 

Having  obtained  A*i>r,  the  first  row  (a*^*)  of  Hc“l  matrix  can  be  obtained 
from  1&-11)  using  c,  fast  Fourier  transform  (FFT)  algorithm  [11]  which  is 
computationally  efficient. 

Now  for  the  elements  of  pth  row  of  He*1  matr'x. 


Since 


g-i2mr  ( i-1)  /N 


(B-13) 


it  can  be  seen  that 


A*  J  2  irr  (p— 1)  /N 
M  l,r 


( B— 14) 


a*P,i  "  a*l,( i-p+1) 


(B-15) 


Furthermore,  it  can  be  derived  from  (&-11)  that 
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a*l,-q  -  a*i>(N_q)  (&-16) 

The  equations  (B-15)  and  (B-16)  indicate  that  Hc"l  is  also  a  circular  matrix, 
and  hence  the  computation  of  its  first  row  is  sufficient  to  determine  all  the 
elements  of  Hc~l. 

Lastly,  it  can  be  easily  verified  that 

a2p-l,2j-l  ■  a2p,2j  -  Real (a*pj- )  (B-17) 

-a2p-l,2j  -  a2p ,2j-l  -  Imag(a*pj)  (B-18) 


where 

aP,j  “ 

In  summary,  Hc“^-  is  obtained  by  computing  the  steady  state  response  of  the 
equivalent  single  degree  of  freedom  model  for  each  engine  order  excitation.  If 
each  blade  sector  is  represented  by  a  greater  number  of  degrees  of  freedom; 
e.g.,  three  in  [6],  the  dimensions  of  matrix  Hc  will  be  3N  x  3N.  In  this  case, 
the  inversion  of  matrix  Hc  will  involve  the  computation  of  the  steady  state 
response  of  the  equivalent  three  degrees  of  freedom  system  for  each  engine  order 
excitation.  In  other  words,  N  number  of  3  «  3  matrices  have  to  be  inverted. 

This  approach  results  in  a  significant  reduction  in  computing  time  with  respect 
to  that  required  for  direct  inversion  of  a  3M  x  3N  matrix,  particularly  when  N 
is  large. 
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Appendix  C 

Computation  of  E(up)  and  E(up2) 

In  the  previous  paper  [4],  formulae  for  computing  E(up)  and  E(up2)  have 
been  developed  when  the  distributions  of  6k j  are  Gaussian.  Here,  these  formulae 
are  generalized  to  include  the  cases  of  uniform  and  triangular  distributions  of 
6k  j. 

First,  the  following  results  are  derived  for  triangular  (Figure  2)  and 
uniform  distributions  (Figure  3): 


E(<5kj) 


2An  6n/2 


fcr  odd  n 


-  '[h+1)  (n+2T*  f°r  even  n  and  triangular  distribution 


(C-l) 


An  3n/2 

;  for  even  n  and  uniform  distributions 


where 

E(<5k2)  -  A2  (C— 2) 

Second,  it  is  recalled  that  the  distribution  of  6k j  is  independent  of  that 
6ks(j#s).  Therefore, 


E(5k*  6kP)  -  E(6kf)  E(6kP)  (C-3) 

Analogous  to  expressions  in  [4],  E(up)  is  computed  on  the  basis  of  first  three 
terms  in  (13)  and  E(UpUq)  only  on  the  basis  of  terms  up  to  6k j ^  (or  6k j 2  6k £2). 
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Hence,  using  (C-l)  and  (03) ,  the  following  expressions  are  obtained  for 
computing  E(up)  and  E(UpUq): 


2  N- 

E(u_)  -  A*  1  w, 
p  j-i 


P.j.j 


(C— 4) 


e(upV  ’a2  j,  VA.j 


N  N 


+  [0  l  yq  ..  i  .-t  •  +  I  l  t  ,(y  •  .  ,+  y  •  .  .+  y  .  .  .}]A 


N  N 


+  j-i  yp»jJ.jtqJ+  J1tq»i(yp,j,j,A+  yP,j.i,j+  yp,i.j.j^A 


N  N 


+  L0  j »iwp »J . jwq . j » j +  I^1wp»j»jwq»i»i+  wp»j»^wq,j,i+wq,i,j^A 

JM  (C-5) 


where 

3.0  for  Gaussian  6k j 
8  *  2.4  for  triangular  ok j 
1.8  for  uniform  6k,- 


(C-6) 


It  is  obvious  the':  E(up^)  can  also  be  computed  from  (C-5)  by  substituting 


q-p- 


TABLE  1 .  System  Parameters  (SI  units) 


=»  0.0114 
kt  -  430000. 
kC  -  45430. 
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SECTION  2 


CALCULATING  THE  STATISTICS  OF  THE  MAXIMUM  AMPLITUDE 
OF  A  MISTUNED  BLADED  DISK  ASSEMBLY 


Nomenclature 


Symbol  - 

Ai 

A, 

c 

cdf 

Cov(y.z) 

diag 

E(y) 

exp 

fo 

f(yi . yn) 

H 

5H 

i 

I 

kC 

kj 

kt 

mt 

Max(Aj/At) 

N 

pdf 

Pr(  ) 

t 


Meaning 

amplitude  of  xj  (mistuned  system) 
resonant  amplitude  of  the  tuned  system 
modal  damping  (tuned  system) 
cumulative  distrioution  function 
covariance  of  random  variables,  y  and  z 
diagonal  elements  of  a  diagonal  matrix 
expected  value  of  random  variable,  y 
exponentia, 

magnitude  of  excitation  force 
joint  pdf  of  random  variables  yi  ,y2,  ... .  yn 
tuned  system  matrix  defined  in  equation  (9) 
matrix  defined  in  equation  (10) 

V=T 

identity  matrix 

structural  coupling  stiffness  between  adjacent  blades 

modal  stiffness  of  jth  blade  (mistuned  system) 

modal  stiffness  of  each  blade  (tuned  system) 

modal  mass  of  each  blade  (tuned  system) 

Max(Ai,A2,...,AN)/At,  maximum  value 

number  of  blades 

probability  density  function 

probability 

time 

jth  element  of  vector  U 
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ul 

i 

Var(y) 

xj 

Pr 

5kj 

Yy,z 

(0 

¥j 

0y 

c 


Uj  for  the  tuned  system  (8kj  =  0) 

variance  of  random  variable,  y 

response  of  jth  blade  (mistuned  system) 

interblade  phase  angle  of  rth  engine  order  excitation 

=  kj  -kt 

correlation  coefficient  between  y  and  z  random  variables 
mean  of  a  random  variable  y 
excitation  frequency  (rad/sec) 
phase  angle  of  Xj 

phase  angle  of  excitation  force  on  the  jth  blade 
standard  deviation  of  random  variable  y 
damping  ratio 
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1.  Introduction 

Mistuning  in  a  biaded  disk  assembly  refers  to  the  variation  in  the  modal  properties  of 
blades  because  of  manufacturing  tolerances.  It  has  received  a  wide  attention  in  the  existing 
literature  [1-3]  because  even  a  small  amount  of  mistuning  can  lead  to  a  large  variation  of 
blades'  amplitudes  within  the  same  assembly.  In  fact,  the  amplitude  of  vibration  of  the  worst 
blade  may  often  be  as  high  as  one  and  one-half  times  the  amplitude  of  a  blade  of  a  perfectly  tuned 
system.  Consequently,  the  fatigue  life  of  a  blade  on  a  mistuned  assembly  may  be  significantly 
lower  than  that  predicted  on  the  basis  of  a  perfectly  tuned  system.  Furthermore,  the  deviations 
of  the  modal  properties  of  blades  from  their  nominal  values  are  statistically  distributed.  As  a 
result,  there  exists  a  large  number  of  possible  sets  of  blades'  modal  parameters  that  a  biaded 
disk  assembly  can  have.  Since  the  maximum  amplitude  will  be  different  for  each  such  set  of 
blades'  modal  parameters,  it  is  important  for  the  designer  to  know  the  statistics  of  the  forced 
response,  particularly  the  probability  that  the  maximum  amplitude  on  any  disk  would  exceed  a 
critical  value  in  an  efficient  manner. 

For  a  given  distribution  of  modal  parameters,  the  statistics  of  the  blades'  amplitudes  can 
be  computed  by  the  following  two  approaches:  1)  numerical  simulations  and  2)  analytical 
technique.  For  numerical  simulations  [4,5],  the  modal  parameters  of  each  blade  of  a  mistuned 
assembly  are  chosen  randomly  from  a  distribution  with  specified  values  of  mean  and  variance. 
Since  the  steady  state  response  of  the  resulting  system  can  be  easily  obtained,  the  statistics  of 
blades'  amplitudes  can  be  generated  by  considering  a  large  number  of  rotor  stages  with  different 
sets  of  blades'  modal  parameters.  However,  this  approach  is  cumbersome  and  expensive  in 
terms  of  computer  time. 

The  analytical  approach  [6-13]  yields  the  statistics  of  the  forced  response  of  a  mistuned 
system  without  performing  numerical  simulations  and,  consequently,  is  computationally 
efficient.  The  technique  developed  by  Huang  [6]  can  yield  the  mean  and  the  variance  of  a  blade’s 
amplitude.  However,  this  analysis  is  valid  only  for  rotor  stages  with  closely  spaced  blades.  In 
addition,  there  was  no  attempt  to  characterize  the  probability  density  function  (pdf)  of  a 
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blade's  amplitude.  Recently,  Sinha  and  Chen  [7-1 0]  have  developed  a  technique  to  compute  the 

statistics  of  the  forced  response  of  a  mistuned  bladed  disk  assembly..  This  technique  yields  the 
pdf  of  a  blade's  amplitude  {Aj)  through  the  joint  statistics  of  Aj  cos<j>j  and  Aj  sinpj  where  <t>j  is 

the  phase  of  the  jth  blade.  This  method  is  valid  for  Gaussian  and  non  Gaussian  distributions  of 
modal  parameters  as  well.  However,  the  distribution  of  the  maximum  amplitude  on  a  bladed 
disk  assembly  was  riot  computed. 

In  this  paper,  an  alternative  analytical  technique  is  developed  to  directly  obtain  the  pdf 
of  the  amplitude  of  a  blade's  vibration  without  using  the  join'  statists  of  Aj  cosfy  and  Aj  sin<})j. 
The  amplitude  of  each  blade  is  expanded  as  a  Taylor’s  series  in  terms  ot  the  Deviations  of  the 
modal  parameters.  Then,  it  is  shown  using  the  central  limit  theorem  that  the  distribution  of 
each  amplitude  is  Gaussian  irrespective  of  the  nature  of  modal  parameter  distributions.  On  the 
basis  of  this  result,  an  efficient  technique  has  been  developed  to  compute  the  distribution  of  the 
maximum  amplitude  on  a  bladed  disk  assembly.  The  comparison  with  the  results  from 
numerical  simulations  indicates  that  the  analytical  technique  is  accurate. 

First,  the  development  of  the  analytical  technique  is  presented.  The  model  for  the  bladed 
disk  assembly  is  shown  in  Figure  1  and  is  exactly  the  same  as  used  in  previous  papers  [7-10]. 
Next,  the  results  of  the  analytical  technique  are  compared  with  those  from  numerical 
simulations. 

2.  Governing  Equation  and  Steady  State  Response 

In  the  model  of  the  bladed  disk  assembly  (Fiugre  1),  the  vibration  each  blade  is 
represented  by  a  single  mode,  mt  and  kj  are  the  modal  mass  and  stiffness,  respectively.  The 
mechanical  coupling  between  adjacent  blades  due  to  the  disk  flexibility  is  represented  by  a 
spring  with  stiffness  kC.  The  viscous  damping  coefficient  c  is  associated  with  each  blade  to 
represent  the  energy  dissipation  from  structural  and  aerodynamic  sources. 

The  mistuning  of  each  blade  is  modeled  by  the  deviation  in  the  blade's  stiffness  alone.  The 
external  force  on  each  blade  is  assumed  to  represent  a  particular  engine  order  and,  as  a  result, 
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is  sinusoidal  in  time  and  differs  in  phase  by  a  constant  amount  from  blade  to  blade  [11].  Under 
the  rth  engine  order  excitation,  the  differential  equation  of  motion  for  the  jth  blade  can  be 
described  as 

+  cxj  +  k.\j  +  kC(x , -  x^j)  +  kCCXj -  xj+1)  =  f0cos (cot  +  yj  ( 1  ) 

where  Y3  =  (j“  l)Pr.P  r  =  27tr/N.:  Because  of  the  periodicity  of  the  rotor  stage  shown  in  Figure 
1,  (j-1)  is  taken  to  be  N  if  j=1,  and  (j+1)  is  taken  to  be  1  if  j=N.  Note  that  kj  =  kt  +  8kj 
where  Skj  is  a  random  variable  with  zero  mean..  Also,  the  distribution  of  8kp  is  independent  of 
8ki  (j*p)..  There  is  no  assumption  made  on  the  nature  of  the  distributions  of  Skj. 

The  steady  state  response  of  the  mistuned  system  under  sinusoidal  excitation  can  be 
described  by  the  harmonic  motion  x(: 

AjCC'ifon  +  tig  (2) 

Substituting  (2)  into  (1)  and  then  equating  the  coefficients  of  coscot  and  siriwt  on  both 
sides,  a  system  of  2N  algebraic  equations  can  be  obtained  as  follows: 

(h,  +  5kj)Ajcos<j)j  +  h3Ajsin  ^  -  h2(A  j+icos  0  J+I  +  AH .cos  r^,,)  =  f0cos  y}  (3) 

(hj  +  SkpA^in  (}>  i  -  hjA^os  (Jr  -  h2(Aj+1sin  <[>j+1  +  A^sin  f0sin  y.  ( 4  ) 

where 

hj=  kt  +  2kC -mtco2  ,  h2  =  kC  ,  h3  =  - cco  (5) 

Rewriting  (3)  and  (4)  in  matrix  form, 

[H  +•  5H]U  =  Fa  (6) 

where 

T 

U  =  (A1cos<fr,,A1sin  (^..^.cos^  ANsin  4>N)  (7) 


Fa  =f0(cos\|/,,sin  y?  ...cosy  y  sin  \|/  N)T 


(8) 


39 


H  = 


f  hj 

h3 

-h2 

0 

0 

0 

...  -  h2 

0 

"h3 

h  i 

0 

-h2 

0 

0 

...  0 

-h2 

-K 

0 

h  3 

-h 

.  0 

....  0 

0 

0 

h  2 

“h3 

h, 

0 

-h2 

.....  0 

0 

"h2 

0 

0 

0 

0 

0 

....  h, 

h3 

1° 

~h2 

0 

0 

0 

0 

....  -h3 

h, 

O) 


The  matrix  H  corresponds  to  the  tuned  system  and  is  circular  in  nature  [7],  The  matrix  5H 

contains  the  deviations  in  the  modal  stiffness  due  to  mistuning  and  is  defined  as  follows:. 

8H  =  diag  (8kp  8k,,  8k2,8k2,..,  8kN,  8kN)  ( q) 

3.  Development  of  the  Analytical  Technique 

3.1  Computation  of  pAj.  oAj.  yAp,  Aq,  and  the  cdf  of  Aj 

The  amplitude  of  a  blade  is  expressed  as  a  perturbation  from  the  tuned  system  (kj=kt, 

Aj=At)  case  via  Taylor  series.  Truncating  the  fourth  and  higher  order  terms: 

n  dA.  ,  N  n  d2A. 

Aj  =  A;  +  X  ^lT8ki  +  2lZ  Z  d  k  d  k  ,5k.5ki 


»=i 


N  N  N 


d’A 


+  -&LL  dk.dk  dk„  8ki8k^8k 


(11) 


i=l  ^=lm=l 


where  5kj  =  kj-kt 


For  any  give  model,  partial  derivatives 


dA.  d2  A. 

i  _ L 


d  k. '  dk.d  k. 


d3  A. 

“d3Mk)k.  “V 


can  be 


computed  in  a  recursive  manner  (Appendix  A).  Having  obtained  these  partial  derivatives,  the 
mean  (pAjh  the  variance  (aAj).  and  correlation  coefficient  between  Ap  and  Aq  (yAp,Aq)  can  be 
directly  computed  (Appendix  B).  From  equation  (11),  Aj  is  a  combination  of  a  large  number  of 
random  variables  which  are  not  highly  dependent  [8].  Therefore,  the  distribution  of  Aj  can  be 
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approximated  as  Gaussian  for  a  large  number  of  blades  (N)  by  the  central  limit  theorem  [12]; 
i.e. 


1 


f(A)  —  /r — 

J  o4v  2n 


exp 


-(A^hAj)  K2a‘A) 


(12) 


n 

Any  linear  combination  of  amplitudes;  i.e.,  ^  aj^j  where  a,  are  arbitrary  numbers, 

j=i  J 


will  also  be  a  combination  of  a  large  number  of  random  variables  which  are  not  highly 
dependent.  Consequently,  it  can  also  be  approximated  to  be  Gaussian.  This  implies  [13]  that  the 

joint  distributions  f(A-j,A2 . An)  of  amplitudes  is  Gaussian;  i.e. 

f(A  A  ,An)  =  -t-'- ---exp-B"  'B 

■J (2lt)  hi  (13) 

where  A  =  [Ai  ,A2 . An] 

B  =  A  -  E(A) 

(Tl)ij  =  ^(Aj  Aj)  -  E(Aj)E(Aj) 

Iti  I  =  det  -q 

Furthermore,  it  should  be  noted  that  the  distributions  of  Ajcosty  and  Ajsinty  have  also  been 

found  to  be  Gaussian  [8];  and  Aj  depends  on  AjCOS<t>j  and  AjSinty  in  a  nonlinear  fashion  as  follows; 

/  2  2  \U2 

Aj  =  l  AjCos2^.  +  AjSin2^.]  (14) 

2 

However,  it  has  been  found  that  E(Aj)  and  E(Aj)  can  be  predicted  with  a  good  accuracy  by 
linearizing  equation  (14)  around  the  tuned  system  response;  i.e., 

Aj  “  Al  +  ’2*A1~I,J(AJC0S<I,J  “  U2j-l) 

+ lA;,'!(Aisini|i-u2,)  ns) 

The  expression  (15)  indicates  that  Aj  can  be  approximated  as  a  linear  combination  of  Ajcos<t>j 

and  Ajsin<j>j  for  the  typical  level  of  mistuning.  Hence,  the  result  that  the  distributions  of  Aj, 
AjCos<)>j  and  AjSin^  are  Gaussian  can  be  explained. 
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3.2  Distribution  of  the  Maximum  Amplitude 

Since  the  failure  of  a  rotor  stage  often  results  from  one  broken  blade  which  corresponds 
to  the  maximum  amplitude,  the  most  important  information  for  the  designer  is  to  know  the 
statistics  of  the  maximum  amplitude  on  a  disk.  The  cdf  of  the  maximum  amplitude  among  the 
given  N  blades'  amplitudes  (Aj)  on  a  disk  can  be  expressed  as  follows 

Pr(Max(A1,A!,..,AN)  £MA) 

=  Pr(A,£MA,A!SMA,..,AN<MJ 

M 

Ma  Ma  * 

“  JQ  10  •••  J  f(Aj,A2,.., An)<9  Ajd  A  2  ...  d  An,  ( -|  e 

o 

where  MAis  a  constant. 

As  discussed  in  section  3.1,  the  joint  pdf  f(Ai,  A2 . An)  can  be  approximated  to  be 

jointly  Gaussian.  Since  E(Aj)  and  the  covariance  matrix  q  can  be  analytically  computed 
(Appendix  B),  one  can  theoretically  evaluate  the  integral  in  equation  (16).  However,  in 
practice,  the  number  of  blades  on  a  disk  is  often  large,  e.g.  N  =  56  [11],  N  =  72  [4-5]. 
Therefore,  a  hyper-dimensional  integration  like  equation  (16)  for  a  general  function  f(Ai,  A2, 
.... .  An)  will  become  extremely  expensive  or  may  be  impossible  in  terms  of  the  computational 
cost.  Hence,  in  order  to  efficiently  evaluate  the  integral  in  (16),  further  assumptions  are 
imposed  on  the  pdf,  f(Ai,  A2 . An). 

If  the  effects  of  the  off-diagonal  terms  of  the  covariance  matrix,  q,  can  be  neglected,  the 

N-dimensional  integration  (16)  can  be  reduced  to  N  one-dimensional  integrals  as  follows 

N  ma 

Pr(Max (Aj»A  2,..,A  N)  <  MA)  =  nj  Af(A)dA.  (17 

j=t  0 

Furthermore,  the  pdf  of  each  Aj  is  identical  as  each  kj  has  an  identical  distribution  [14]. 
Therefore,  equation  (17)  can  be  further  simplified  as  follows 

Pr(Max(A,,A r..,A  „)  <  MA>  =  (j“*f(A,)dA ,)  ( ,  g 

Comparing  equations  (16)  and  (18),  it  can  be  seen  that  an  almost  impossible  N- 
dimensional  integration  has  been  reduced  to  a  simple  one-dimensional  integration.  The  validity 
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of  equation  (18)  is  established  in  the  next  section.  It  should  be  noted  that  the  distribution  of  the 
maximum  amplitude  is  not  Gaussian. 

4.  Comparison  with  Numerical  Simulations 

In  order  to  verify  the  accuracy  of  the  analytical  technique,  the  analytical  results  are 
compared  with  those  from  numerical  simulations.  For  numerical  simulations,  the  steady  state 
responses,  AjCos4>j  and  AjSin<j>j,  of  the  mistuned  bladed  disk  assemblies  are  obtained  by  solving 
the  2N  simultaneous  linear  equations  (3)  and  (4).  Then,  the  statistics  of  Aj  and  Max(Aj)  are 
computed  from  a  large  sample  (arbitrarily  chosen  to  be  3000)  of  bladed  disk  assemblies  with 
different  sets  of  deviations  in  blades'  modal  parameters,  8kj.  The  values  of  8kj  are  chosen  from 
a  Gaussian1  distribution  with  specified  means  and  variances.  The  IMSL  [15]  subroutine 
DRNNOA,  has  been  used  to  generate  the  Gaussian  distributions. 

The  number  of  blades  is  chosen  to  be  twenty-four.  The  excitation  frequency  is  chosen  to 
be  6328.8  rad/sec,  which  corresponds  to  the  resonant  condition  for  the  tuned  system  under  the 
third  engine  order  excitation  (r=3).  The  third  engine  order  corresponds  to  the  45  degree 
interblade  phase  angle  of  the  excitation  forces.  The  nominal  values  of  system  parameters  are 
shown  in  Table  1  and  they  are  representative  of  a  high  pressure  turbine  stage  [4].  The  damping 
value,  c  =  1 .443  Newton-sec/meter,  corresponds  to  about  1 .0%  damping  factor,  which  is 
typical  of  those  observed  in  gas  turbine  engine  rotor  stages  from  aerodynamic  and  structural 
sources  [5]. 

In  order  to  test  the  accuracy  of  the  analytical  results,  the  standard  deviations  (crskj)  of 

the  mistuned  stiffnesses  are  chosen  to  be  7500  Newton/meter  and  10000  Newton/meter, 
respectively.  Also,  E(8kj)=0. 


1  The  distributions  of  amplitudes  have  been  found  to  be  insensitive  to  the  types  of  modal 
parameters'  distributions  [9]. 
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4.1  Distribution  of  Amplitudes 

From  Figure  2,  it  can  be  seen  that  »he  analytical  results  compare  well  with  those  from 
the  numerical  simulation,  even  though  the  accuracy  of  the  analytical  result  decreases  as  the 
value  of  CTgkj  increases. 

4.2  Distribution  of  the  Maximum  Amplitude 

As  shown  in  Figure  3,  the  distribution  of  the  normalized  maximum  amplitude 
Max(Aj/At)  (At  is  the  resonant  amplitude  of  the  tuned  system)  from  the  analytical  technique, 
equation  (18),  compares  well  with  that  from  the  numerical  simulation  for  ogkj  =  7500. 
However,  tor  askj  -  1 0000,  the  accuracy  becomes  smaller.  This  is  due  to  the  error  in 
estimating  the  distribution  of  the  amplitude  by  the  analytical  technique,  equat'on  (12),  for  ogkj 
=  10000.  It  shoul;  ,oe  noted  that  the  analytical  approximation,  equation  (18),  is  found  to  be 
quite  accurate  if  the  distribution  of  the  amplitude  from  the  numerical  simulation  is  used,  Figure 
3(b). 

The  magnitudes  of  the  off-diagonal  terms  in  the  covariance  matrix,  t\  can  be  evaluated  by 
examining  the  correlation  coefficients,  YAp.Aq-  Because  of  the  periodicity,  it  has  been  found 
that  the  values  of  YAm.Aj  depend  on  I  m-jl .  Hence,  only  the  values  of  yai  >Aj  are  presented  here. 
As  shown  in  Figure  4  for  each  askj,  the  correlation  coefficients  (yAi.Aj)  oscillate  around  zero 

and  decrease  as  the  distance  between  two  blades  increases.  Also,  the  absolute  values  of  most  of 
yAi  ,Aj  are  small  and  almost  half  of  them  are  very  close  to  zero.  Even  though  some  of  yAi.Aj  are 
not  zero,  the  cumulative  effects  of  the  off-diagonal  terms  in  the  covariance  matrix,  q,  are  small 
and  the  accuracy  of  the  analytical  technique,  equation  (18),  is  acceptable  (Figure  3). 

As  mentioned  earlier,  the  correlation  coefficients,  YAi.Aj,  decrease  from  the  nearer 
neighbours  to  the  farther  ones,  Figure  4.  This  implies  that  the  magnitudes  of  the  off-diagonal 
terms,  YAm.Aj,  become  smaller  as  I  m-jl  becomes  larger.  Hence,  one  should  expect  the  accuracy 
of  the  analytical  results  to  be  better  for  the  rotor  stage  with  a  larger  number  of  blades,  while 
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worse  for  that  with  a  smaller  number  of  blades.  To  study  this  effect,  the  number  of  blades  is 
chosen  to  be  N  =  8  and  N  =  56.  The  standard  deviation  of  8kj  is  equal  to  7500. 

As  shown  in  Figure  5,  relatively,  most  of  the  yAi.Aj  are  not  zero  for  N  =  8,  while  for  N 
=  56  more  than  four-fifths  of  the  ya-|  , Aj  are  close  to  zero.  Hence,  the  accuracy  of  the 
distribution  of  the  Max(Aj/At)  from  the  analytical  technique  is  worse  of  N  =  8,  while  much 
better  for  N  =  56,  Figure  6.  Since  the  number  of  blades  is  often  large  in  practice  (e.g.  N  =  56 
[11],  N  =  72  [4,5],  the  analytical  approximation,  equation  (18),  may  be  successfully  applied 
in  these  practical  situations. 

5  Conclusions 

To  efficiently  compute  the  statistics  of  the  maximum  amplitude  of  a  mistuned  bladed  disk 
assembly,  an  analytical  technique  has  been  developed  in  this  paper.  The  analytical  technique 
uses  the  direct  Taylor  series  expansion  in  terms  of  the  perturbations  in  an  amplitude  as  a 
function  of  perturbations  in  modal  stiffnesses.  This  technique  can  compute  the  mean,  the 
variance,  the  cdf  of  an  amplitude,  and  the  correlation  coefficients  between  any  two  amplitudes. 
Also,  the  cdf  of  the  maximum  amplitude  on  a  rotor  stage  can  be  computed. 

The  validity  vf  the  analytical  technique  has  been  corroborated  by  comparison  with  the 
results  from  the  numer  cal  simulation.  Also,  the  analytical  technique  developed  in  this  paper  is 
^alid  for  Gaussian  and  non-Gaussian  mistuned  parameters  as  well.  The  CPU  time  for  the 
analytical  technique  is  found  to  be  much  smaller  than  that  for  numerical  simulations.  On  a  VAX 
8550  computer,  the  CPU  time  for  numerical  simulations  and  the  analytical  technique  are  1079 
seconds  and  9.59  seconds,  respectively  for  N  =  56  and  c§kj  =  10000. 
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APPENDIX  A 

d  A.  d 2  A  d\ . 

COMPUTATION  OF  _1  U-  AND  ^ L — 

Tk~  ’  d  k  dk*  dkdk^dk, 


From  (7),  the  amplitude,  Ai,  can  be  represented  as 


Ar(“l,  +  “2/ 


(A.  1 ) 


Continually  differentiating  (A.1), 


dAi  if  9  %-i  d*2l') 

TiT  =  7^u%-»“3TEr  +  u^Tir  J  *  (A. 2) 

a2A  gAj^A,  d  U2j-1  a  u2._,  du2j  d  u2j 

<9  k  <9  k  A. I  d  k  d  k.  +  d  k  d  k  +  d  k  ’  d  k 

t  *  A  x  i  x  >  x  i 

9  U2j-1  9  U2j 

+  u2i-i<9k,dk  +u2j-^kTF  (A-3> 

i  >  4  1/ 


and  for  p  >  2, 


s,</  s2  ...  d  sp 


■3P-1 

5  S. 


UzH  d  Sj<?  s2 ...  <9  sp  u2j^  s,<9  s2  ...  <?  s 


<H  “q+1 
"I 


(A. 4) 


where  si,S2 . sp  represent  a  combination  of  p  variables  chosen  from  N  variables,  kj.  And, 


91=  Aj,g2  =  U2j-1,  -  U2j,sign(gi )  =  -1,  sign(g2)  =  sign(g3)  =  1.  From  equation  (6), 

d\J  (H  +  8H)  cT’u 


Sld  S2  d  Sp 


--(H  +  8H)_1X 


+  5H) _ d^U _ 

sq  dsl...dsrlds^...dsp' 


(A. 5) 


where  3(H+8H)/3Sq  can  be  obtained  for  a  given  model  of  the  bladed  disk  assembly.  For  the 


model  shown  in  Figure  1, 


d  (H  +  5H) 


diagj’O...  I'  1...0J 


(A. 6) 


=  aP,q> 


where  sq  =  kj.  If  (H”1)Piq 
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u  u  2m -1  _ 

(9  k.  U2j-ia2m-Uj-l  “  U2ja2m-Uj 

<9  u2„, 

(9  k  ~  ~  U2j-ia2mij-l  “  U2ja2m^j 


(A. 7) 


(A. 8) 


dk.dk, 

J  2 


£?  k  (9  k  ,  [ a  2m -1,22-1  (9  k.  +a2m-1.22^k 

J  2  V  J  J 

9  U  2j-l  3  U2^ 

+  a2m-1.2j-l  ^  k,  +  a2">-1.2j  (9  k,  j 

Sa„~,  , ,  in„  ,  £  %-■  ,  £V| 

2m. .22-1  <?k  a2m.2.<<9k  &  2”  .2i-‘  <9  k  .  +  &  2m-2J  <9  k  .  j 

J  J  "  *  / 


(A. 9) 


(A. 10) 


dV. 


^2u2m 


<9k<9k<9k.  I  a2m-ui-i  (9  k  <9  k  +  a2m-ui^  k  d  k.  +  a2m-i,2j-i  d  k  <9  k 

2  J  >  V  2  J  2  J  2  > 

a*-*  Aa  ihu  ) 

+  a2m-l.2id  k,(9  k  +ai«-UM9k.9k  +  a2m-U2  <9  k  <9  k  /A  11\ 

2i  j  i  j  i  /  ( A .  1  1  ) 

^U2m  f  ^«2i-l  .  d^2i  . 


(9  k  (9  k  <9  k.  ~  I  a2mii-i^  k  9k.  a2m.2id  k  (9  k.  a  2<n.2J-i  d  k  d  k 

2  J  1  V  2  J  2  J  2  1 


+  a2m,2j^k,^k  +a2m.21-l^  k  <9  k.  +a2n.^l^  k  k  }  (A-12) 

2i  j  i  J  i  / 

It  should  be  noted  that  all  the  partial  derivatives  are  evaluated  at  the  tuned  system  (kj  = 


kt)  .esponse.  Furthermore,  since 


<9%  <92A, 

<9  k.<9  k  “  d  k.d  k.’ 

i  2  2  i 


(A. 13) 


d3  A.  <93A.  <93A 

_ i _ j _ > 

d  k.d  k  <9  k  ~  d  k.d  k_<9  k  “  d  k  d  k.<9  km 

i  2  ra  i  m  2  2  i  m 

<93  Aj  <93A.  <93A; 

”  <9  ky9  km<9  ki  ~  d  k  md  k.d  k  x~  d  k  md  k  xd  k.’ 

only  one  term  has  to  be  computed  in  the  above  two  groups,  (A.13)  and  (A.14). 


(A. 14) 
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APPENDIX  B 

CALCULATION  OF  A  >  AND  y. 

Ai  i  aj,a> 


From  equation  (11), 


n  A  ,  N  N 

a  =  a  +y  4^6k  +  sk 

p  1  ^  d  k  1  2!  9  k  d  k  1  > 

•  i=i  j=i  i  j 


1=1 


N  N  N 


iD2 


d'A, 


3!^f-  <9  k  <9  k  ^  k:Sk*6ki6k^* 

1=1  j=l  i=l  1  J  A 


Then  the  above  equation  can  be  expressed  as 

Ap  -  At  =  tp  +  Wp  +  yp 


where 


•p'X'Jsk,) 

j=i 

__  d  Ap 
1pj 

J 

wp=ZZ  wP.j*<5kj><5kk) 

j=lk=l 


w 


1 


d\, 


pj-k  2!  <9k.dk/ 


yP  =  XIX  yPJ,k.i(5kJ)(5kk)(5kP 

j=l  k=l  J=l 


1 


<93Ar 


p.jJt  3!  <9  k  ,<9  k.<9  k 

l  k  j 


Then  HA  ,  c\  and  yA  ^  can  be  computed  from  Ref.  [8]. 


j  > 


(21) 

(B.1) 


(B.2) 


(B.3) 


(B  .4) 
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Table  1  Data  for  Simulation  (SI  Unit) 


mt  =  0.0114  kC  =  45430.0 
kj  =  430000.0  c  =  1.443 
f0  =  1.0 
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Numerical  Simulation 
Analytical  Technique  [8] 
Analytical  Technique 


Figure  2  cdf  of  Aj:  Numerical  Simulation  and  Analytical  Technique  N  *=  24,  c 
45°,  (a)  ask]  =  7500  (b)  051. j  =  10000 


=  1 .443,  (Jr 


Cumulative  Probability 
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Max(Aj / At) 


Max(A:-/At) 


Figure  3  cdf  of  Max{A]/A{):  Numerical  Simulation  and  Analytical  Technique  N  =  24,  c  = 
1.443,  pr  -  45°,  (a)  agkj  -  7500  (b)  cr5kj  «  10000 


Blade  Number 


Figure  5*  Correlation  Coefficients}  f 

c  =  1.443,  fir  =  45°,  o6k,  =  7500  6kf.  Gaussian 

(a)  N  =8  (b )N  =56 


Numerical  Simulation 
Analytical  Technique 


Figure  6.  cdf  of  Max(Aj/At)  .  ,  ,  . 

Numerical  Simulation  and  Analytical  Technique 
c  =  1.443,  fiT  =  45°,  cskj  =  7500  6kf.  Gaussian 

(a)  N  =8  (b )N  -56 


SECTION  3 

ANALYTICAL  STATISTICS  OF  FORCED  RESPONSE  OF  A  MISTUNED 
BLADED  DISK  ASSEMBLY  IN  SUBSONIC  FLOW 
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Abstract 

The  statistics  of  the  forced  response  of  a  structurally  and  aerodynamically  coupled  bladed 
disk  assembly  have  been  computed  efficiently  by  the  analytical  technique.  The  validity  of  the 
analytical  technique  has  been  corroborated  by  comparison  with  the  results  from  numerical 
simulations.  Lastly,  the  effects  of  the  following  parameters  on  the  statistics  of  the  maximum 
amplitude  have  been  studied:  aerodynamic  coupling,  fluid  density,  and  stagger  angle. 


Nomenclature 

■Symbol 

h 

A‘ 

J 

At 

b 

c 

cdf 

diag 

i 

kC 

kj 

kt 

■^whr.^hhr 

mt 

N 

S 

U 

V 


Meaning 

amplitude  of  Xj  (mistuned  system) 

=  A^xp  * 1 

(resonant)  amplitude  of  xtj  (tuned  system) 
semicord  length 

viscous  damping  (tuned  system) 
cumulative  distribution  function 
diagonal  elements  of  a  diagonal  matrix 

structural  coupling  stiffnesses  between  adjacent  b'ades 
modal  stiffness  of  jlh  blade  (mistuned  system) 

modal  stiffness  of  each  blade  (tuned  system) 

nondimensional  aerodynamic,  coefficients 

modal  mass  of  each  blade  (tuned  system) 

number  of  blades 

spacing  between  adjacent  blades 

vector  defined  in  equation  (10) 

velocity  of  fluid  flow 

response  of  jth  blade  (mistuned  system) 
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response  of  j,h  blade  (tuned  system) 
stagger  angle 

interblade  phase  angle  of  r,h  engine  order 
=  kj  -  kf 

location  of  elastic  axis 
excitation  f-equency  (rad/sec) 
phase  angle  of  Xj 
fluid  density 

standard  deviation  of  random  variable  y 
damping  ratio 
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1.  Introduction 

Mistuning  in  a  bladed  disk  assembly  refers  to  the  variation  in  the  modal  properties  of 
blades  because  of  manufacturing  tolerances.  It  has  received  wide  attention  in  the  existing 
literature  [1-5]  because  even  a  small  amount  of  mistuning  can  lead  to  a  large  variation  of 
blades'  amplitudes  within  the  same  assembly.  In  fact,  the  amplitude  of  vibration  of  the  worst 
blade  may  often  be  as  high  as  one  and  one-half  times  the  amplitude  of  a  blade  of  a  perfectly 
tuned  system.  Consequently,  the  fatigue  life  of  a  blade  on  a  mistuned  assembly  may  be 
significantly  lower  than  that  predicted  on  the  basis  of  a  perfectly  tuned  system.  Furthermore, 
the  deviations  of  the  modal  properties  of  blades  from  their  nominal  values  are  statistically 
distributed.  As  a  result,  there  exists  a  large  number  of  possible  sets  of  blades’  modal 
parameters  that  a  bladed  disk  assembly  can  have.  Since  the  maximum  amplitude  will  be 
different  for  each  such  set  of  blades'  modal  parameters,  it  is  important  for  the  designer  to  know 
the  statistics  of  the  forced  response,  particularly  the  probability  that  the  maximum  amplitude 
on  any  disk  would  exceed  a  critical  value  in  an  efficient  manner. 

To  study  the  dynamic  response  of  a  mistuned  bladed  disk  assembly,  the  aerodynamic 
effects  have  been  often  simuated  as  constant  viscous  damping  [6-11]  associated  with  each 
blade's  motion.  This  simple  representation  of  aerodynamic  forces  provides  a  good  basis  for  the 
development  of  physical  concepts  and  mathematical  techniques  which  are  required  to  sludy  the 
effects  of  mistuning  on  the  forced  response  of  a  bladed  disk  assembly.-  However,  this  model  does 
not  consider  certain  important  features  of  cascade  aerodynamics,  e.g.,  the  aerodynamic  coupling 
between  each  blade  and  the  dependence  of  aerodynamic  loads  on  the  iriterbiade  phase  angle,  the 
frequency  of  blade  motion,  the  Mach  number,  and  various  cascade  parameters.  The  effect  of 
limiting  aerodynamic  and  structural  coupling  has  been  studied  by  Basu  and  Griffin  [12],  and 
the  statistics  of  amplitudes  were  obtained  by  the  numerical  simulation.  The  aerodynamic 
coupling  between  each  blade  results  in  the  dense  system  equations  which  require  a  large  amount 
of  computation  to  solve  them.  As  a  result,  to  generate  the  statistics  of  amplitudes  by  the 
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numerical  simulation,  the  computational  cost  for  an  aerodynamically  coupled  system  is  much 
higher  than  that  of  a  structurally  coupled  system.  Consequently,  the  first  objective  here  is  to 
efficiently  calculate  the  statistics  of  amplitudes  of  a  mistuned  assembly  whose  bladed  are 
coupled  structurally  and  aerodynamically.  Then,  the  next  objective  is  to  investigate  the  effects 
of  various  cascade  parameters  on  the  statistics  of  the  amplitudes  of  a  mistuned  assembly. 

The  spring/mass  model,  Figure  1,  is  used  in  this  paper  to  represent  a  rotor  stage 
[6-11].  This  model  considers  only  one  mode  of  vibration  for  each  blade.  The  unsteady 
aerodynamic  loads  due  to  motions  of  blades  in  a  rotor  stage  have  been  calculated  by  the 
linearized  cascade  theory  [13-15]  for  the  subsonic  flow.  For  computing  aerodynamic  loads, 
the  rotor  stage  is  modeled  as  a  two-dimensional  cascade  of  blades,  Figure  2.  The  blade  is 
represented  by  a  "typical  section"  motion  which  is  assumed  to  vibrate  in  the  bending  mode, 
Figure  3, 

In  addition  to  the  interblade  phase  angle  (Pr  =  2rcr/N),  the  aerodynamic  forces  also 

depend  on  the  stagger  angle  (a),  the  ratio  of  blade  spacing  to  cord  length  (S/2b),  the  reduced 
velocity  (cob/V),  the  Mach  number,  and  the  location  of  the  elastic  axis  (u). 


2.  Governing  Equations  and  Steady  State  Response 

As  stated  in  the  Ref.  [3]  each  blade  has  identical  amplitude  and  the  interblade  phase  angle 
is  conslant  in  each  mode  for  the  tuned  cascade.  Hence,  for  the  rtfl  tuned  mode,  the  response  of 

iL 

;he  jm  blade,  x{j,  can  be  represented  as 


xtj=  A' exp 


i(jp  ,+toi) 


(1) 


where  \  jS  a  complex  number.  Its  magnitude  represents  the  amplitude  of  motion.  In  a 

mistuned  assembly  for  which  each  blade  has  a  different  response,  the  motion  of  the  jth  blade, 
Xj,  can  be  represented  as  a  liner  combination  of  all  the  tuned  modes  corresponding  to  all  the 


interbiade  phase  angles  [3];  i.e. 


N-l 


N-S 


r -0  rU) 


(2) 
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where  hor  are  complex  coefficients.  Thus,  for  the  mistuned  assembly,  the  equation  of  motion 
«r 

for  the  jth  blade  is  [3] 

mlxj  +  cx,+  kjXj  +  kCCXj-  x;_,)  +  kC^ -  xj+1) 

=  Lrcpb  [i  hhr-y-  +  *  WJ  exp  ,(“>t+jP')  (  3  ) 

r=0 

where  j  =  1,2,  ...  ,  N.  L  and  p  are  unit  length  and  fluid  density,  respectivelv.  The  stiffnesses 
kj(=kt  +  5kj)  are  chosen  as  random  variables  to  repiesent  the  mistuning  in  this  study.  The 
.Jhhr  is  the  nondimensional  lift  coefficient  [3]  for  the  rtfl  tuned  mode  and  ^whr  is  the  rth 

nondimensionai  lift  coefficient  due  to  the  wakes  shed  from  upstream  periodic  obstructions. 
These  nondimensional  coefficients  (ihhr.^whr)  are  obtained  from  the  aerodynamic  code,  NASA 
Lewis  Research  Center  [15].  The  linear  viscous  coefficient,  c,  represents  the  damping  from 
the  structural  sources.  The  value  of  the  viscous  damping  coefficient,  c,  is  selected  so  that  the 
total  damping  in  the  tuned  model  is  one  percent  of  the  system’s  critical  damping  which  is 
usually  found  in  practice  [12]. 

For  the  steady  state  response,  the  motion  of  the  j,h  blade  of  *he  mistuned  system  can  be 
represented  as  a  harmonic  motion 

Xj=  exp‘“'  (4) 

C  1* 

where  A;(  =  A^xp  J)  is  the  complex  amplitude.  Comparing  equations  (2)  and  (4),  one  can 


obtain 


N-l 


a'=X  mV 


2itrj/N) 


r=0 


(5) 


Then  using  the  inverse  discrete  Fourier  transformation  [8],  equation  (5)  yields 

J=1 


(6) 


Now  substituting  equations  (4)  and  (6)  into  (3), 

(-  mico2  +  icci)+  kt  +  2kC  +  8kJ  A'  -  kCA'^  -kCA 


j+i 


L  rcpb  o) 

FT 


2„2  f  N  TN-l 

I  V 


I  I  ^hhreXPl2ra<J",+1)/N|A‘[  =  Lftpb3©2]^  ^whiexp’ 

.  «=1  Lr=0 


N-l 


i2nrj/N 


(7) 


J  ) 


The  above  equation  can  be  written  in  tlv»  matrix  form  as  follows 
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[Ha  +  5K]U=Fa  (8) 

where  the  matrix  HA  and  the  vector  FA  are  defined  according  to  equation  (7).  And, 

T 

U  =  (A,  cos 0 ,,  AjSin  ^)p..,AN  cosrj)^  ANsin  <|)N)  (  9  J 

6H  =  diag(8k1,  5k1(  8k2,  5k2,.., ... Sk^  8kw),  (10) 


In  order  to  understand  the  nature  of  the  mistuned  response,  the  wake  excitation  (,ewhr) 

in  one  particular  engine  order  will  be  considered  at  a  time  throughout  this  study..  As  one  can 
see  in  equation  (7),  the  motion  of  each  blade  is  coupled  to  all  the  blades  by  the  aerodynamical 
terms  which  result  in  the  dense  system  equation.  Consequently,  the  computation  time  for  the 
numerical  simulation  gets  enhanced. 

For  a  tuned  system  (8k.- 0  and  A'=  A°mexp  l2,CT(j"m)/N  [3])  under  the  rth  engine 
order  excitation,  equation  (7)  can  be  simplified  as  follows: 

[-  m(co2  +  k,  +  2kC(l  -  cos(?.7r:/N)  -  L7tpb2to2i hhr  +  icoc]  A'  =  LrcpfrVi  whr  exp l2,trj/N, 

(11' 

The  resonant  frequency  (cor)  and  the  damping  ratio  (Q  of  the  tuned  system  can  be  computed  as 
follows: 

'kt  +  2kC(l  -  cos(2jct/N))1y 

0)  r  =  5  »  ii? 

L  mt  +  L7cpb2Real(ihhr)  J  U  ^ 

c-7ipb2corImag(ihJ 

~  2(mi  +  L7cpb2Real(ihhr)o)r)  ‘  (13 

3.  Analytical  Technique  to  Compute  thr*  Statistics  of  Blades'  Amplitudes 

Even  though  the  system  with  aerodynamic  coupling  is  more  complicated  than  that  with 
only  structural  coupling,  it  still  remains  linear.  Hence,  the  analytical  technique  developed  in 
Ref  [11]  can  be  directly  applied  to  cor  pute  the  statistics  of  blades'  amplitudes:  the  mean,  the 
variance  and  the  cdf  of  an  amplitude,  the  correlation  coefficients  between  any  two  amplitudes, 
and  the  cdf  of  the  maximum  amplitude  on  a  rotor  stage.  Therefore,  tk  objective  here  is  to 
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study  the  accuracy  of  this  technique  and  the  aerodynamic  effects  0,1  the  statists  of  tb'> 
maximum  amplitude, 

4,  Comparison  with  the  Numerical  Simulation 

To  verify  the  applicability  of  the  analytical  techniques  to  the  mistuned  system,;  with 
structural  and  aerodynamical  coupling  between  blades,  the  cdf  of  Max(Aj/A0  from  the 
analytical  technique  [11]  is  compared  with  that  from  the  numerical  simulation,  where  At  is 

the  resonant  amplitude  of  the  tuned  system.  For  the  numerical  simulation,  the  IMSL 
subroutine  DRNNOA  [16]  is  used  to  choose  the  modal  stiffnesses  (5kj)  from  the  Gaussian 

population  with  E[5k;]=0  and  csi<j=7500  Newton/meter.  The  response  of  the  resulting  bladed 
disk  assembly  is  obtained  by  !  'ing  the  N  linear  algebraic  equations  (7)  in  terms  of  Ajexp1^ 
with  the  IMSL  subroutine  DL2LCG  [16].  The  cdf  of  Max(Aj/At)  of  each  blade  is  then  generated 
by  considering  a  large  number  (arbitrarily  chosen  to  be  3000)  of  bladed  disk  assemblies  with 
different  sets  of  blades’  modal  stiffnesses.  The  nominal  values  of  the  system  parameters  are 
shown  in  Table  1,  where  the  aerodynamic  parameters  are  chosen  from  Ref.  [3].  The  engine 
order  of  excitation  is  chosen  to  be  seven;  i.e.  r  =  7,  The  amount  of  damping  from  the 
aerodynamic  source  in  the  tuned  model  is  found  to  be  0.73  percent  of  the  system's  critical 
damping  under  the  seventh  engine  order  excitation.  Hence,  the  vaiue  of  the  viscous  damping 
coefficient,  c,  has  been  chosen  to  be  0.22  percent  of  the  critical  damping;  i.e.  c  =  0.3193 
Newton-s  ;c/meter.  Since  the  semicord  length,  b,  on  the  right  hand  side  of  equation  (7)  acts  as 
a  scaling  factor,  it  is  chosen  to  be  unity  in  this  study. 

In  Figure  4,  the  cdf  of  Max(Aj/At)  from  the  analytical  technique  [3]  agrees  well  with 
that  from  the  numerical  simulation.  This  result  indicates  that  the  analytical  technique 
developed  in  Ref  [11]  can  also  efficiently  and  accurately  compute  the  statistics  of  blades’ 
amplitudes  of  a  mistuned  assembly  whose  blades  are  coupled  structurally  and  aerodynamically. 
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5  Paramedics  Studies 

5.1  The  Effect  of  Aerodynamic  Couplings  Among  Blades 

The  statistical  response  of  the  system  developed  in  Ref.  [11]  whose  blades  are  only 
couplet?  structurally  with  the  immediate  neighbours  is  compared  with  that  of  the  sysiem 

developed  in  this  paper,  in  each  system,  the  equivalent  viscous  damping  ration  (Q  is  equal  to 
1.0%  The  cumulative  disiributon  functions  of  Max(Aj/Aj)  indicate  that  the  aerodynamical^' 

coupled  system  has  a  higher  probability  of  having  a  larger  amplitude.  Figure  E.  These  results 
indicate  that  the  system  in  which  the  aerodynamic  effects  are  modeled  as  <i  vise;  us  damping 
associated  with  each  c.sde's  motion  undei  estimates  the  effects  of  mistiming. 

5.2  The  Effect  of  Fluid  Density 

From  equation  (7).  it  is  cisa*  that  the  fluid  density  centrals  the  aerodynamic  damping 
and  the  coupling  between  blades.  The  statistics  of  amplitudes  are  computed  tm  the  Allowing 
values  of  the  fluid  censity:  OJSp,  o,  i.2Gp  and  l.op.  ihe  cdf  of  Max(Aj/At)  for  each  value  of 
the  fluid  density  is  shown  in  Figure  6.  It  indicates  that  the  system  with  a  larger  value  of  the 
fluid  density  has  a  smaller  probability  of  having  a  higher  amplitude.  In  other  words,  the 
impact  of  mistuning  is  more  severe  on  the  system  with  a  lower  value  of  the  fluid  density. 

5.3  The  Effect  of  Stagger  Angle 

The  aerodynamic  coefficients  {-inhr  and  iWhr)  depend  cn  the  stagger  angle  of  the  cascade 
(Figure  2).  As  a  result,  the  aerodynamic  damping  and  the  coupling  between  blades  will  be 
different  for  each  stagger  angle.  The  statistics  of  amplitudes  are  computed  for  the  following 
values  of  the  stagger  angle  (a):  15°,  30*,  45°,  60°  and  75s.  Evert  though  the  variation  in  the 
equivalen*  iscous  damping  ration  {£)  is  large  for  these  stagger  angles  (£=0.80%,  0.83%, 
0.91%,  1.08%,  1.54%  for  a=15°,  30°,  45°,  60®  and  75°„  rospectively),  the  results,  cdf  of 
Max(Aj/At),  indicate  the  the  statistics  of  the  blades'  amplitudes  are  not  highly  sensitive  to  the 
variation  in  the  stagger  angle,  Figure  7. 
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From  the  above  results,  it  is  clear  that  the  damping  ratio  alone,  which  oniy  depends  on 
jihhr.  equation  (13).  is  not  enough  to  explain  the  mistuning  effect.  The  effect  of  other 
aerodynamic  coupling  terms  (ih.hj.  j  *  r)  should  be  considered  as  well.  To  consider  the  impact 
cf  these  terms,  one  can  compute  the  inverse  va  ue  of  the  condition  number  (t)  of  the  system 
matrix.  Ha  The  hverse  value  of  the  condition  number  can  measure  how  close  a  matrix  is  to 
me  singularity  [17].  The  smaller  the  value  of  r  is,  the  closer  is  the  matrix  to  the  singularity, 
and  the  higher  order  terms  will  be  consequently  more  important.  The  values  of  t  for  the 
system  matrix  with  stagger  angle  <x=15°,  30°,  45°,  60°  and  75°,  are  0.04356,  0.04385, 
0.04471,  0.04676,  and  0.0445,  respectively.  The  results  indicate  that  the  degree  with 
which  each  system  is  close  to  the  singularity  is  almost  the  same.  Consequently,  the  statistical 
distribution  of  amplitudes  are  identical  in  spite  of  the  variation  in  the  equivalent  viscous 
damping  for  the  rlfl  tuned  mode. 

In  section  5.1,  t=0.05878  for  the  structurally  coupled  system  [11],  while  t=*0. 04596 
tor  the  aerodynamically  coupled  system.  Therefore,  even  though  the  total  damping  ratio  is  1  % 
for  the  rtfl  tuned  mode  in  each  case,  the  scatter  in  amplitudes  of  the  aerodynamically  coupled 
system  is  greater  than  that  of  the  structurally  coupled  system  (Figure  5). 

In  section  5.2,  -1=0.0383,  0.04596,  0.05319  and  0.06  for  the  fluid  density  equal  to 
0.5  p,  p,  1.25p,  and  1.5p,  respectively.  And,  the  corresponding  equivalent  viscous  damping 
ratios  are  0.80%,  1.00%,  1.20%  and  1.39%,  respectively.  Therefore,  the  system  matrix  is 
less  singular  and  the  equivalent  viscous  damping  ratio  also  increases  as  the  fluid  density 
increases.  These  tacts  explain  why  the  system  is  less  sensitive  to  mistuning  as  the  fluid 
density  increases  (Figure  6). 

6.  Conclusions 

The  statistics  of  the  forced  response  have  been  computed  by  the  analytical  technique  [11] 
fcr  a  linear  assembly  whose  blades  are  coupled  structurally  and  aerodynamically.  The  validity 
of  the  analytical  technique  has  been  corroborated  by  comparison  with  the  results  from  the 
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numerical  simulation.  The  CPU  time  for  the  analytical  technique  is  found  to  be  much  smaller 
than  tnat  for  the  numerical  simulation.  On  a  VAX  8550  computer,  the  CPU  times  for  the 
numerical  simulation  a  id  the  analytical  technique  are  3068.0  seconds  and  9.72  seconds, 
respectively. 

The  parametric  studies  indicate  the  following: 

1  )  The  mistuning  effect  on  a  blade's  amplitude  may  be  underestimated  by  the 

structurally  coupled  system  in  which  the  aerodynamic  effects  are  modeled  as  a 

viscous  damping  associated  with  each  blade. 

2  )  The  probability  of  having  a  blade  with  a  larger  amplitude  decreases  as  the  fluid 

density  increases. 

3  )  The  stagger  angle  of  the  cascade  has  little  impact  on  the  statistics  of  a  blade's 

amplitude  under  the  subsonic  flow  condition. 
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Table  1:  Data  for  Simulation  --  Aerodynamically  Coupled  System  (SI  Unit) 


mt/(Lnpb2)  -  258.5 
cob/V  =  0.55 
v  =  0.5 
mt  =  0.0114 
kt  =  430000 
c  =  0.3193 
E[5kj]  =  0 


M  =  0.7 
a  =  54° 

S/(2b)  =  0.534 
kC  =  45430.0 
N  =  56 

tTSkj  =  7500.0 
6kj:  Gaussian 


:  viscous  damping  from  structural  source  only 


Figure  1  .  Spring/Mass  Model  of  a  Rotor  Stage  (Linear  System) 


Figure  2,  Cascade  of  Blades 


Equivalent  Section  of  a  Blade 


Figure  3, 


Cumulative  Probability 


Numerical  Simulation 
Analytical  Technique 


4.  cdf  of  Max(Aj /At)  . 

Numerical  Simulation  and  Analytical  Technique 
N  -  56,  c  -  0.3193,  pr  =  45°,  a6kj  -  7500,  6k}-:  Gaussian 
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Figure  6. 


The  Effect  of  Fluid  Density 
on  the  cdf  of  \'I ax{Aj  <  A:) 

N  =  56,  c  =  0,3193,  (3r  “  45° , 


o(k}  =  7500,  6kj\  Gaussian 


Figure  7: 


The  Effect  of  Stagger  Angle  , 

on  the  cdf  of  Max\AdAt J  Q  Gaussian 

N  —  56,  c  =  0.3193,  /?r  =  45°,  *«*,-  “  <oUU’ 


